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Abstract 

This research report presents an extension of Cumulative of Choco constraint solver, 
which is useful to encode over-constrained cumulative problems. This new global con- 
straint uses sweep and task interval violation-based algorithms. 

1 Introduction 

The Cumulative global constraint provides a pruning algorithm which takes account of all 
activities at the same time, which has been proved to be much more efficient than CP 
approaches considering a conjunction of more primitive constraints. This representation 
as a global constraint has been widely studied in the literature and integrated into many 
constraint systems [TJ [2 [3l |H [5] . 

Definition 1 Let A = {a\, . . . ,a n } denote a set of n non-preemptive activities. For each 
a G A, 

• start[a] is the variable representing its starting point in time. 

• dur[a] is the variable representing its duration. 

• end[a] G end is the variable representing its ending point in time. 

• res[a] is the variable representing the discrete amount of resource consumed by a, also 
denoted the height of a. 

Definition 2 Consider one resource with a limit of capacity max-capa and a set A of n 
activities. At each point in time t, the cumulated height ht of the set of activities overlapping 
t is h t = J2aeAMart[a]<t<endla] min{D(res[a})) . 

The Cumulative global constraint Jljj enforces that: 

• CI: For each activity a £ A, start [a] + dwr[a] = end [a]. 

• C2: At each point in time t, h t < max_capa. 

In this research report, we deal with cumulative over-constrained problems that may 
require to be relaxed w.r.t. the resource capacity at some points in time, to obtain solu- 
tions. This motivates the definition of a new global constraint, dedicated to over-constrained 
instances of problems. 
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2 Cumulative Constraint with Over-Loads 



This section presents the Sof tCumulativeSum constraint, useful to express our case-study 
and deal with significant instances. D(x) = {min(D(x)), . . . ,max(D(x))} denotes the do- 
main of variable x. 

2.1 Pruning Independent from Relaxation 

The Sof tCumulativeSum constraint that will be presented in section [2~2l implicitly defines a 
Cumulative constraint with a capacity equal to max_capa. To prune variables in start, sev- 
eral existing algorithms for Cumulative can be used. We recall the two filtering techniques 
currently implemented in the Cumulative of Choco solver [6], which we extend to handle 
violations in section l2~2l 

2.1.1 Sweep Algorithm [4] 

The goal is to reduce domains of start variables according to the cumulated profile, which is 
built from compulsory parts of activities, in order not to exceed max-capa. This is done in 
Choco using a sweep algorithm [4]. The Compulsory Part [7] of an activity is the minimum 
resource that will be consumed by this activity, whatever the final value of its start. It is 
equal to the intersection of the two particular schedules that are the activity scheduled at 
its latest start and the activity scheduled at its earliest start. As domains of variables get 
more and more restricted, the compulsory part increases, until it becomes the fixed activity. 

Definition 3 The Compulsory Part CP (a) of an activity a is not empty if and only 
if max{D(start\a])) < min(D(end[a])). If so, its height is equal to min(D(res[a])) on 
[max(D(start[a])), min(D{end[a]))\ and null elsewhere. 




Activity A 

non empty compulsory part 



Activity B 

empty compulsory part 



i i i i i i i i i i . m 



Figure 1: Compulsory parts. 

Example 1 Let a be an activity such that start[a] € [1,4], with a duration fixed to 5, and b 
be an activity such that start[b] G [1,6], with a duration fixed to 3. As depicted by Figured 
activity a has a non empty compulsory part, while activity b has an empty compulsory part. 

Definition 4 The Cumulated Profile CumP is the minimum cumulated resource consump- 
tion, over time, of all the activities. For a given t, the height of CumP at time t is equal 
to 

D(min(res[a])) 

a£A/t£[max(D(start[a])),min(D(end[a]))[ 

That is, the sum of the contributions of all compulsory parts that overlap t. 
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Next figure shows an example of a cumulative profile CumP where at each point in time 
t the height of CumP at t does not exceed maxjeapa. 




max_capa 



Figure 2: Cumulated profile. 



Algorithm The sweep algorithm moves a vertical line A on the time axis from one event 
to the next event. In one sweep, it builds the cumulated profile incrementally and prunes 
activities according to this profile, in order not to exceed maxjzapa. An event corresponds 
either to the start or the end of a compulsory part, or to the release date of an activity. All 
events are initially computed and sorted in increasing order according to their date. Position 
of A is 5. At each step of the algorithm, a list ActToPrune contains the activities to prune. 

• Compulsory part events are used for building CumP. All events at date 5 are used 
to update the height sumh of the current rectangle in CumP, by adding the height if 
it is the start of a compulsory part or removing the height otherwise. The first event 
with a date strictly greater than 6 gives the end 6' of the current rectangle in CumP, 
finally denoted by ([6, S'[, sumh}- 

• Events corresponding to release dates d such that 6 < d < 5' add some new activities 
to prune, according to ([6,S'[, surrih) and max-capa (those activities which overlap 
([5, S'[, sumh)). They are added to the list ActToPrune. 

For each a £ ActToPrune that has no compulsory part registered in the rectangle 
([5, 5'[, sumh), if its height is greater than max-capa — sunih then the algorithm prunes its 
start times so this activity doesn't overlap the current rectangle of CumP. Then, if the due 
date of a is less than or equal to 5' then a is removed from ActToPrune. After pruning 
activities, S is updated to 5' . 

Pruning Rule 1 Let a e ActToPrune, which has no compulsory part recorded within the 
rectangle ([5, S'[, sumh) ■ If sumh+min(D(res[a])) > max-capa then ]5 — min{D(dur\a])), S'[ 
can be removed from D(start[a\). 

Time complexity of the sweep technique is 0(n * log(n)). Please refer to the paper for 
more details w.r.t. this algorithm [4]. 

2.1.2 Energy reasoning on Task Intervals [3 12 [9] 

The principle of the energy reasoning is to compare the resource necessarily required by a set 
of activities within a given interval of points in time with the available resource within this 
interval. Relevant intervals are obtained from starts and ends of activities ("task intervals"). 
This section presents the rules implemented in Choco. 

Notation 1 Given ai and dj two activities (possibly the same) s.t. min(D(start[ai])) < 
max(D(end[aj])), we denote by: 
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• I(a,i,a ) the interval [min(D(start[ai])),max(D(end[aj]))[. 

• S(ai,a ) the set of activities whose time-windows intersect I( ai . aj ) o-nd such that their 
earliest start is in I( ai . aj ), that is, S( a4i0 .) = {a G A s.t. min(D(start[ai})) < 
min(D(start[a])) < max(D(end[aj]))} . 

• AREA( 0ii0 .) the number of free resource units in /( aj)a .), that is, AREA( aii0 .) = 
(max(D(end[aj])) — min(D(start[ai]))) * max-capa. 

Definition 5 A lower bound W( ai;a .)(a) of the number of resource units consumed by a G 
s (ai, aj ) on I( aiiaj ) is 

min(D(res[a])) 
* 

min[min(D(dur[a\)),max(0,max(D(end[aj])) — max(D(start[a])))] 
Feasibility Rule 1 IfY^aes, , w (a i ,a j )( a ') > AREA (tt ) then fail. 

Algorithm The principle is to browse, by increasing due date, activities aj G A 
and for a given a,j to browse, by decreasing release date, activities a, G A such that 
min(.D(s£ari[aj])) < m&x(D(end[aj])). Hence, at each new choice of a, or aj more ac- 
tivities are considered. For each couple (a,i,a,j), the algorithm applies the feasibility rule 
ffl 

Suppose the activities sorted by increasing release date i.e. mm(D(start[ai])) < 
mm.(D(start[a,2])) < ■ • • < mhi(D(start[a n ])), then: 

• If min(D(start[ai-i])) = mm(D(start[ai])) then 

^(oi-i.oj) = ^(ai,^) and therefore Z)oes (o , ^W^a^Ca) = 

• Else 

We have S( 0i _ 1)0 ,) = S( 0i , 0j ) U {a fe G < i-1 A min(£>(start[a fe ])) = 

min(D(sfart[aj_i]))} i.e. we add all activities with same release date than 
activity Oi_i. Hence, E a6S(0i _ i aj) ^.^^(W = Eaes^., ^a^O*) + 

S{fe<j-lAmin(£)(start[o fc ]))=min(r>(stori[a i _ 1 ]))} ^(oijOj ) ( a ) • 

Therefore, the complexity for handling all intervals Ii ai ,a ) IS 0(n 2 ). 
2.2 Pruning Related to Relaxation 

Specific constraints on over- loads exist in real- world applications. To express them, it is 
mandatory to discretize time, while keeping a reasonable time complexity for pruning. These 
specific constraints are externalized because they are ad-hoc to each application. On the 
other hand, the following constraints capture the generic core of this class of problems. 

1 By definition, W( a . }a .){a) is independent of a< so W( a ._ lia .)(a) = W( a ., a .) (a). 
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• Sof tCumulative extends Cumulative of Choco by maintaining over-loads variables at 
each point in time, and by pruning activities according to maximal available capacities 
given by upper bounds of these variables instead of simply considering max-capa. This 
constraint can be used with various objective criteria. 

• The Sof tCumulativeSum extends the Sof tCumulative. It is defined to deal more 
efficiently with a particular criterion: minimize the sum of over-loads. 

2.2.1 SoftCumulative Constraint 

Definition 6 Let A be a set of activities scheduled between time and m, each consuming a 
positive amount of the resource. SoftCumulative augments Cumulative with a second limit 
of resource ideal jcapa < max-capa, and for each point in time t < m an integer variable 
costVar[t\. It enforces: 

• CI and C2 (see Definition [2p. 

• C3: For each point in time t, costV ar[t] = max{Q, h t — ideal jcapd) 

costVar[t]: maximum values of domains 



max_capa 













i 











ideal_capa 



costVar[i] lower bounds: 2 2 2 

Figure 3: Example of a SoftCumulative constraint. 

Example 2 Figured depicts an example of a cumulative profile CumP where at each point 
in time t the height of CumP at t does not exceed the maximum value in the domain of 
its corresponding violation variable costVar[t\. Time points \, 2, 3 are such that CumP 
exceeds ideal jcapa by two. Therefore, for each point, the minimum value of the domain of 
the correspondong variable in costVar should be updated to value 2. 

Next paragraph details how the classical sweep procedure can be adapted to the 
SoftCumulative constraint. 

Revised Sweep pruning procedure. The limit of resource is not max-capa as in the 
Cumulative constraint. It is mandatory to take into account upper-bounds of variables in 
costVar. One may integrate reductions on upper bounds within the profile, as new fixed 
activities. However, our discretization of time can be very costly with this method : the 
number of events may grow significantly. The profile would not be computed only from 
activities but also from each point in time. 
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We propose a relaxed version where for each rectangle we consider the maximum costVar 
upper bound. This entails less pruning but the complexity is amortized: the number of time 
points checked to obtain maximum costVar upper bounds for the whole profile is m, by 
exploiting the sort on events into the sweep procedure. 

Pruning rule [2] reduce domains of start variables from the current maximum allowed 
over-load in a rectangle^ 

Pruning Rule 2 Let a G ActToPrune, which has no compulsory part recorded within the 
current rectangle. If sumh+min(D(res\a\y) > ideal_capa+max t ^!s^[( max {D(costVar[i\))) 
then ]S — min(D(dur[a])), S'[ can be removed from D(start[a\). 

Proof: For any activity which overlaps the rectangle, the maximum capacity is upper-bounded 
by ideal-capa + max(D(costVar[t]),t £ [S, S'[). 

Time complexity is 0(n * log(n) + m), where m is the maximum due date. 

Revised task Intervals energy reasoning. This paragraph describes the extension of 
the principle of section 12. 1.21 to the Sof tCumulative global constraint. 

Feasibility Rule 2 AREA( ajia .) = 

ideal jzapa + max(D(costVar[t])) 

tE [min(D (start [a;])) ,max(D(end[aj])) [ 

J /Eaes, , Wf ai>a ) {a) > AREA (oij0 ) then fail. 

Proof: At each time point t there is ideal_capa + max(D(costV ar[t])) available units. 

Efficiency of rule [T] can be improved by this new computation of AREA( a . aj ) (the previ- 
ous value was an over-estimation). Since activities ai are considered by decreasing release 
date, it is possible to compute incrementally AREA( a . a .y Each upper bound of a variable 
in costVar are considered once for each aj. Time complexity is 0(n 2 + n * m), where m is 
the maximum due date. 

Next paragraph explains how minimum values of domains of variables in costVar are 
updated during the search process. 

Update costVar lower-bounds. Update of costVar lower bounds can be directly per- 
formed within the sweep algorithm, while the profile is computed. 

Pruning Rule 3 Consider the current rectangle in the profile, [S,5'[. For each t G [5, S'[, 
if suvrth — idealjcapa > min(D(costV ar[t])) then [min(D (costVar [t])), sumh — ideal_capa[ 
can be removed from D(costVar[t]). 

Proof: From Definitions [3] and 0] 



2 The upper bound max(D(costVar[t])) is the maximum value in the domain D(costVar[t]). Since these 
variables may be involved in several constraints, especially side constraints, the maximum value of a domain 
can be reduced during the search process. 
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Usually the constraint should not be associated with a search heuristic that forces to 
assign to a given variable in costVar a value which is greater than the current lower bound 
of its domain. Indeed, such a search strategy would consist of imposing at this point in 
time a violation although solutions with lower over-loads at this point in time exist (or even 
solutions with no over- load). However, it is required to take into account this eventuality 
and to ensure that our constraint is valid with any search heuristic. If a greater value is 
fixed to a variable in costVar, until more than a very few number of unfixed activities exist, 
few deductions can be made in terms of pruning and they may be costly (for a quite useless 
feature). Therefore, we implemented a check procedure that fails when all start variables 
arc fixed and one variable in costVar is higher than the current profile at this point in time. 
This guarantees that ground solutions will satisfy the constraint in any case, with a constant 
time complexity. 

2.2.2 SoftCumulativeSum Constraint 

Definition 7 SoftCumulativeSum augments SoftCumulative with an integer variable 
cost. It enforces: 

• CI and C2 (see Definition [2]), and C3 (see Definition [Sp. 

• The following constraint: cost — X)te{o m-i} costV ar[f\ 

Pruning procedures and consistency checks of SoftCumulative remain valid for 
SoftCumulativeSum. Additionally, we aim at dealing with the sum constraint efficiently 
by exploiting the semantics. We compute lower bounds of the sum expressed by cost vari- 
able. Classical back-propagation of this variable can be additionally performed as if the sum 
constraint was defined separately. 

Example 3 The term back-propagation is used to recall that propagation of events is not 
only performed from the decision variables to the objective variable but also from the objective 
variable to decision variables. For instance, let x\,X2 and x% be 3 variables with the same 
domain: Vi G {1, 2, 3}, D{xi) — {1,2,3}. Let sum be a variable, D(sum) = {3,4, ...,9} 7 
and the following constraint sum = 2 3} x i- Assume that 2 is removed from all D(xi). 

The usual propagation removes values 4, 6 and 8 from D(sum). Assume now that all values 
greater than or equal to 5 are removed from D(sum). Back-propagation removes value 3 
from all D(xi). 

Sweep based global lower bound. Within our global constraint, a lower bound for the 
cost variable is directly given by summing the lower bounds of all variables in costVar, 
which are obtained by pruning rule [3l These minimum values of domains were computed 
from compulsory parts, not only from fixed activities. 

LBx = min(D(costVar[t})) 

t£{0,...,m-l} 

LB\ can be computed with no increase in complexity within the sweep algorithm. 

Interval based global lower bound. The quantity X^es Wr aiia .)(a) used in feasi- 
bility rule [T] provides the required energy for activities in the interval I( ai ^ aj )- This quantity 
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may exceed the number of time points in I( ai .aj) multiplied by ideal jzapa. We can im- 
prove LB\, provided we remove from the computation over- loads yet taken into account in 
the cost variable. In our implementation, we first update variables in costVar (by rule [3]), 
and compute LB\ to update cost. In this way, no additional incremental data structure is 
required. 

To obtain the new lower bound we need to compute lower-bounds of cost which are local 
to each interval 7( 0il a,-)- 

Definition8 lbu ai ,a) = Stez min(D{costVar[t])) 

Then, next proposition is related to the free available number of resource units within a 
given interval. 

Proposition 1 The number FreeAreA( q . Qj ) of free resource units in I( ai . aj ) s. t. no 
violation is entailed is (max(D(end[aj])) — min(D(start[ai]))) * idealjcapa. 

Proof: From Definition [Jj 

From Definition [5] and Proposition [TJ FreeArea^. aj j +^i( ai ,aj) is the number of time 
units that can be used without creating any new over- load into the interval /( ai ,aj) compared 
with over-loads yet taken into account in lbi^ a . a .y 

Definition 9 lNC (oi , a ) = Y.aes, , W / (a,.a J )(a)-FREEAREA (aij ) - Z&i(„ i)0 ,) 

lNC( aij0j .) is the difference between the required energy and this quantity. Even if one variable 
in costVar has a current lower bound higher that the value obtained from the profile, the 
increase lNC( ai aj .) is valid (smaller, see Definition [S]) . We are now able to improve LB\. 

LB 2 = LBi + max( auaj ) eA 2(lNC^ ai!aj )) 

Another lower bound can be computed from a partition P of [0,m[ in disjoint intervals 
obtained from pairs of activities (a.i,aj): LB^pj = LB\ + ^ 7 £P lNC( a . a .). Obvi- 

ousy LB 2 < LB^py However, time complexity of the algorithm deduced from rule [H in 
the Sof tCumulative constraint is 0(n 2 + n * m). This complexity should reasonably not 
be increased. Computing LB 2 can be directly performed into this algorithm without any 
increase in complexity^ 

Pruning Rule 4 If LB 2 > min(D(cost)) then [min{D{cost)) 1 LB 2 [ can be removed from 
D(cost). 

Proof: From Definition [7] LBi is a lower bound of cost. Since intervals are disjoint, by Definition 
[9] the quantity L_B(p) is a lower bound of cost. LB2 < L-B(p). Therefore LB2 is a lower bound of 
cost. The pruning rule holds. 

Aggregating local violations. Once the profile is computed, if some activities having a 
null compulsory part cannot be scheduled without creating new over-loads, then LB\ can 
be augmented with the sum of minimum increase of each activity. This idea is inspired from 
generic solving methods for over-constrained CSPs, e.g., Max-CSP [10] . Our experiments 
shown that there is quite often a way to place any activity without creating a new violation. 
This entails a null lower bound. Therefore, we removed that computation from our imple- 
mentation. We inform the reader that we described the procedure in a preliminary technical 
report [IT] . 

3 On the contrary, determining a relevant partition P from the activities would force to use an independent 
algorithm, which can be costly depending on the partition we search for. Finally, we decided to use only 
LB 2 . 
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Instance 


cost value 


Sof tCumulative 

+ external sum 


S o f t Cumul at i v e Sum 


1 





92 (0.07 s) 


92 (0.01 s) 


2 


2 


417 (0.29 s) 


94 (0.01 s) 


3 


10 


> 30 s 


63 (0.01 s) 


4 


2 


1301 (0.59 s) 


194 (0.06 s) 


5 


6 


19506 (13.448 s) 


97 (0.01 s) 


6 





53 (0.00 s) 


53 (0.00 s) 


7 


10 


> 30 s 


90 (0.01 s) 


8 


6 


> 30 s 


152 (0.07 s) 



Table 1: Number of nodes of optimum schedules with n — 9, m — 9, durations between 1 and 4, 
resource consumption between 1 and 3, ideal-capa = 3, max-capa — 7. 



2.2.3 Implementation 

Constraints were implemented to work with non fixed durations and resource consumptions. 
Table Q] compares the two constraints on small problems when the objective is to minimize 
cost. Results show the main importance of LB2 when minimizing cost. 

3 Extension 

The global constraint presented in this research report can be tailored to be suited to some 
other classes of applications. If the time unit is tiny compared with the makespan, e.g., 
one minute in a one-year schedule, the same kind of model may be used by grouping time 
points. For example, each violation variable may correspond to one half-day. Imposing 
a side constraint between two particular minutes into a one-year schedule is generally not 
useful. For this purpose, the Sof tCumulative constraint can be generalized, to be relaxed 
with respect to its number of violation cost variables. 

3.1 RelSoftCumulative constraint 

Notation 2 To define RelSoftCumulative we use the following notations. Given a set of 
activities scheduled between and m, 

• mult £ {1,2,..., m} is a positive integer multiplier of the unit of time. 

• Starting from 0, the number of consecutive discrete intervals of length mult that are 
included in the interval [0, ..m[ is \m/mult \ . J = {0, 1, . . . , \m/mult\ — 1} is the set 
of indexes of such intervals. Hence, to each j £ J corresponds the interval [j *mult, j * 
mult + 1, . . . (j + 1) # mult — 1] . 

Definition 10 Let A be a set of activities scheduled between time and m, each consuming 
a positive amount of the resource. RelSoftCumulative augments Cumulative with 

• A second limit of resource idealjcapa < max-capa, 

• The multiplier mult £{1,2,..., m}, 

• For each j G J an integer variable costVar[j\. 
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It enforces: 

• CI and C2 (see Definition HP . 

• C'4: For each j G J , costV ar[j] = Ylt^*muU lt * max (^i h t — ideal_capa) 

Example 4 We consider a cumulative over-constrained problem with n activities scheduled 
minute by minute over one week. The makespan ism = 2940. Assume that a user formulates 
a side constraint related to the distribution of over-loads of resource among ranges of one hour 
(mult — 60,) in the schedule, for instance "no more than one hour violated each half-day". 
The instance of RelSoftCumulative related to this problem is defined with [2940/60], i.e., 
49 violation variables, J = {0, 1, . . . ,48}. For each range indexed by j G J , the constraint 
C3 is: costVar[j] = sum^^J^ 1 max(0, h t — ideal -capa). The side constraint is then 
simply expressed by cardinality constraints over each half-day, that is, each quadruplet of 
violation variables: {costVar[0\, • • • , costyar[3]}, {costl^ar[4], • • • , costVar[7]}, etc. 

It is possible to reformulate rules of section 12.21 to make them suited 
to RelSoftCumulative. Firstly, rule [2] can be re-written for the constraint 
RelSoftCumulative. 

Pruning Rule 5 Let a G ActToPrune, which has no compulsory part recorded within the 
current rectangle. 

If sumh + min(D(res[a])) > ideaLcapa + max(max(D(costVar[j]),j/[d,S'[nJ'(j) ^ 0)j3 
then ]S — min(D(dur[a])), 5'[ can be removed from D(start[a\). 

Similarly, rule [2] is reformulated as follows: 
Feasibility Rule 3 AreA( q . j0 .) = 

ideal-capa + max (D(costVar[[t/ mult \})) 

t£[min(D(start[ai])),rnax(D(end[aj])) [ 

I fJ2 a es iaita . ) w (ai,a,j){a) > Area^^.) then fail. 

To update minimum values of domains of variables in costVar in RelSoftCumulative, 
we simply have to update for each violation variable, during the sweep, the current sum 
of over-loads of its time points. This may be done only by maintaining one value and one 
index, but for sake of clarity we use the following notation. 

Notation 3 Given a set of ranges in time indexed by J = {0, 1, . . . , \m/muV\ — 1}, 
costarray is an array of integers. All of them are initially set to 0. They are one-to- 
one mapped with elements in J . 

Next rule reformulates rule [3] for RelSoftCumulative. 

Pruning Rule 6 Consider the current rectangle in the profile, [S,5'[. For each t G [5, S'[, 
if sumn — ideal-capa > then: 

1. costarray [[t/mu/tj] <— costaxray[\t /mult\] + sum^ — ideal_capa 

2. if costarray[[i/mu^J] > min(D(costVar[\t/mult\])) then: 

The range [min{D{costVar[\t/mult\])), costarray[[t/mw?tJ][ can be removed from 
D(costVar[[t/mult\}). 

4 The range of index j in J corresponding to time point t is [t/mult]. Hence, the set of j such that 
[S,S'[nJ{j) ^ 0)) is {[8/mult\, . . . , [(8' - l)/mult\}. 
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3.2 RelSoftCumulativeSum constraint 

Definition 11 RelSoftCumulativeSum augments RelSoftCumulative with an integer 
variable cost. It enforces: 

• CI and C2 (see Definition ^j), and C4 (see Definition [TO]) . 

• The following constraint: cost = YljeJ costVar[j] 

Sweep based global lower bound, as for the Sof tCumulativeSum constraint presented 
in section 12.2.21 a lower bound for the cost objective variable is given by summing the 
lower bounds of all variables in costVar, without any increase in complexity in the sweep 
algorithm. 

LBi = min(D (costVar 

Interval based global lower bound. The task interval energetic reasoning presented in 
section [2.2.21 remains the same, except the evaluation of the quantity lb±, which corresponds 
to over-loads expressed by variables in costVar array (sec definition [8]). 

Within the filtering algorithm of the Sof tCumulativeSum constraint, lb\ is the sum of 
min(costVar[t]) over all points int time t within each considered interval, providing that, in 
the implementation, all variables in costVar array arc updated before computing lb\. With 
respect to the RelSoftCumulative constraint, by definition 1111 if mult is greater than or 
equal to two then it is not possible to evaluate at each point in time t the exact over-load 
at t. Under-estimating lb\ would be false because this leads to a over-estimation of task 
interval based lower bounds (see Definition [5]) . 

Therefore, we compute an over-estimation of Ibi, the tightest possible according to the 
definition of the constraint. 

Notation 4 Let I = [a,b] be an interval of points in time included in [0,m[, and the set of 
indexes for intervals in time J — {0, 1, . . . , \m/mui\ — 1}. For each j 6 J, #(J, j) is the 
number of time points in common between the range indexed by j and I: 

#(i,j) = \[a,b]nj(j)\ 

We can now reformulate Definition [8] for the RelSoftCumulative. The idea is to eval- 
uate for each intersecting the current interval I( ai ,a )-i the minimum value between 
min(D(costV ar[j])) and the maximum possible over-load in / n J(j), which is equal from 
Definition [IT] to #(/, j) * (max(D(costV ar[j])) — ideal-capa). 

Definition 12 Let I( ai . aj ) be a task interval, and J = {j e J , #i(-f(o»,a )) > 0}- 

lb\ = min(min(D(costVar[j))), #(/, j) * (max(D(costVar[j])) — ideaLcapa)) 
je.J 

4 Conclusion 

This report proposed several filtering procedures for a global Cumulative constraint which 
is relaxed w.r.t. to its capacity at some points in time. We provided the extension of our 
global constraint for the case where side constraints are related to ranges in time which are 
larger than one time unit. 
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